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1 Introduction 



From the perspective of dynamical systems theory a dynamical system shows reac- 
tion type dynamics if its phase space possesses a bottleneck type geometry. The 
system spends a long time in one phase space region before it finds its way through 
a bottleneck to another phase space regions. In the terminology of chemistry these 
phase space regions correspond to 'reactants' and 'products', respectively, and the 
bottleneck connecting the two is called a 'transition state'. For Hamiltonian systems, 
the phase space bottlenecks result from equilibrium points of a certain stability type, 
namely equilibria for which the matrix associated with the linearization has one pair 
real eigenvalues of opposite sign and / — 1 complex conjugate pairs of imaginary 
eigenvalues, where / is the number of degrees of freedom. Such an equilibrium is 
called a saddle x center x • • • xcenter, and we will refer to them as 'saddle' for short. 
Near a saddle the energy surface (locally) bifurcates from an energy surface with 
two disjoint components (the reactants and products) for energies below the energy 
of the saddle to an energy surface consisting of a single connected component for 
energies above the energy of the saddle. The single component energy surface has 
a wide-narrow- wide geometry, and the system has to pass or 'react' through this 
bottleneck in order to evolve from reactants to products or vice versa. The most 
widely used approach to compute reaction rates in chemistry is to place a dividing 
surface in the bottleneck region and compute the rate from the flux through this 
dividing surface. This forms the basis of Transition State Theory that was developed 
by Eyring, Polanyi and Wigner in the 1930's [Il[2]. For this approach to be useful the 
dividing surface needs to have the property that it is crossed exactly once by all re- 
active trajectories (i.e., trajectories moving from reactants to products or vice versa) 
and not crossed at all by the non-reactive trajectories (i.e., trajectories staying on 
the reactants or products side). The construction of a dividing surface which solves 
this so called recrossing problem has posed a major difficulty in the development of 
transition state theory. In the 1970's Pechukas, Pollak and others |3lll] showed that 
for systems with two degrees of freedom, a recrossing free dividing surface can be 
constructed from a periodic orbit (the Lyapunov orbit associated with the saddle). 
The longstanding problem of how to generalize this idea for systems with more than 
two degrees of freedom has been solved only recently [21 IH] , see also [TJ [H] . Here the 
role of the Lyapunov orbit is taken over by a normally hyperbolic invariant manifold 
(NHIM) [HI HU] which forms the anchor for the construction of a dividing surface in 
the general case. The NHIM also gives a precise meaning of the transition state. It 
can be viewed as the energy surface of an invariant subsystem with n — 1 degrees 
of freedom, which as an unstable 'super molecule' [llj is located between reactants 
and products. The NHIM is not only of central significance for the construction of 
a recrossing free dividing surface but also gives detailed information of the geometry 
and the mechanism of reactions. This arises from the fact that the NHIM has sta- 
ble and unstable manifolds which are of one dimension less than the energy surface. 
They form the separatrices which separate the reactive trajectories and non-reactive 
trajectories in the energy surface. In fact they form tubes which snake through the 
phase space and for a phase space point in the region of reactants to be reactive (i.e.. 
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for a point leading to a reactive trajectory if it is taken as an initial condition for 
Hamilton' s equations) it has to be contained in a certain volume that is enclosed 
by the stable manifold of the NHIM. Similarly only points in the region of prod- 
ucts can be reached by reactive trajectories emanating from the reactants if they are 
contained in a certain volume enclose by the unstable manifold of the NHIM. This 
information is crucial not only for the study of state specific reactivities but also 
for the control of reactions. Most importantly for applications, the dividing surface, 
the NHIM and the local pieces of its stable and unstable manifolds can be explicitly 
constructed from a Poincare-Birkhoff normal form expansion about the saddle equi- 
librium point [U [12] . The normal form gives a symplectic transformation to a new 
set of phase space coordinates in terms of which one can give simple formulas for the 
phase space structures. The inverse of the normal form transformation then allows 
one to construct the phase space structures in the original ('physical') coordinates. 

In this paper we address the question of how the phase space structures can 
be constructed for the relative equilibria of a rotational symmetry reduced A^-body 
system. Here we have in mind that the bodies represent the N atoms that con- 
stitute a given molecule (described in the Born-Oppenheimer approximation) but 
the theory presented also works for other A^-body systems like celestial iV-body sys- 
tems. For molecular reactions, the excitation of rotational degrees of freedom almost 
always plays an important role. The phase space bottleneck are then no longer as- 
sociated with equilibrium points but with larger sets of phase space points which 
however become again equilibria (so called relative equilbibria) if considered in the 
rotational symmetry reduced system. In this context the reduction is important 
for both conceptual and computational reasons. We are particularly interested in a 
reduction which facilitates the Poincare-Birkhoff normal construction of the phase 
space structures around a relative equilibrium in the same fashion as in the case 
of a usual saddle equilibrium. The theory of reduction for A^-body systems is well 
developed [13l HH [151 [El [13 IH]- However, although the structure of the reduced 
space is well-known, the explicit choice of suitable coordinates for the reduced space 
remains as a challenging problem [IS]. If the angular momentum is zero, the re- 
duced space is symplectomorphic to the tangent bundle of the shape or internal 
space (see Sec. [s]). For non- vanishing angular momentum however, this is no longer 
the case [13]. This makes it difficult to apply classical techniques such as Poincare- 
Birkhoff normal form [20], and that point is the motivation of our paper (see also 
[211 [22] ) . There exist a huge literature on the analysis of reduced A/'-body systems 
[231[a[2a|26l[271[28llM[30l|3ll[M[MliaiMl[Ml[3^ theory 
is well-developed. Nevertheless a more canonical way of expressing the equations of 
motion seems to be useful, and we hope that the theory presented in this paper will 
also be of interest for studies of A^-body systems other than the reaction dynamics 
studied here. 

In our approach we mainly follow the account given by Littlejohn and Reinsch 
|17j . Their reduction procedure can be viewed as a generalization of the free rigid 
body reduction scheme. Namely, the Hamiltonian is written as the sum of rotational 
kinetic energy, vibrational kinetic energy and the potential by passing to the body 
coordinates. The only remaining coordinates are the body angular momentum, the 
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shape or internal coordinates and the conjugate momenta of the internal coordinates. 
Using the body angular momentum the equations of motion have a non-canonical 
form which is not so suitable for the local analysis near equilibrium points |4 HI42lH3] . 
Making use of the fact that the reduced space is locally a product of the constant 
body angular momentum sphere and the cotangent bundle of the internal space we 
write the equations of motion in a Hamiltonian form in this paper. We do this by 
choosing some suitable coordinates on the angular momentum sphere. 

This paper is organized as follows. In Sections [2] and |3] we review the phase 
space structures governing reaction dynamics across saddle equilibrium points and 
the reduction of the rotational symmetry of A^-body systems, respectively. The main 
result of this paper is contained in Sec. |4] where we introduce a canonical coordinate 
system on the symmetry reduced space which facilitates Poincare-Birkhoff normal 
form computations near equilibrium points. The theory is illustrated for the limiting 
case of a rigid body and the triatomic example in Sec. [5j Conclusions and an outlook 
are given in SecjHj 

2 Phase space structures governing reactions 
across saddles 

In this section we discuss the phase space structures governing reaction dynamics 
near saddle equilibrium points. Before studying the case of a general (nonlinear) 
Hamiltonian vector field it is useful to first consider the linear case. For the details 
we refer to [HI US] . 

2.1 The linear case 

Consider a linear Hamiltonian vector field with a saddle on the phase space IR,'^ x R-^ 
where / > 2 (we comment on the case / = 1 at the end of this section). Consider 
the quadratic Hamiltonian function 



The corresponding linear Hamiltonian vector field has a saddle equilibrium point at 
the origin, i.e. the matrix associated with the linear vector field has the pair of real 
eigenvalues it A and / — 1 pairs of complex conjugate imaginary eigenvalues iiw^, 
k = 2, . . . , f . We define the constants of motion 




(1) 



k=2 



^i=Pi-qi, h=vl + <ll, k = 2,...,f, 



(2) 



which up to positive prefactors agree with the energies in the individual degrees of 
freedoms. 
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Consider a fixed energy E > 0, where is the energy of the saddle. Rewriting 
the energy equation H2{(l, p) = E in the form 

^pl + j:Yipl + 4)=E + ^q! (3) 
k=2 

one sees that each fixed gi G IR defines a topological (2/ — 2)-dimensional sphere. 
The energy surface 

= {(q, P) G : i^2(q, p) = E} (4) 

thus has the topology of a spherical cylinder R x S^^^^ . The 'radius' of the family 
of (2/ — 2)-dimensional sphere ([S]) becomes smallest for qi = 0, and this in fact can 
be used to define a recrossing free dividing surface (see the introduction). Setting 
= on the energy surface gives the (2/ — 2)-dimensional sphere 

^DS~' = {(q^ P) e : ^2(q, p) = , gi = 0} . (5) 

The dividing surface Sjjg divides the energy surface into the two components which 
have gi < (the 'reactants') and gi > (the 'products'), respectively, and as qi = 
dH2/dpi = \pi 7^ for pi 7^ the dividing surface is everywhere transverse to the 
Hamiltonian flow except for the submanifold where qi = pi = 0. For qi = pi = 0, 
the energy equation ([s]) reduces to ^^{=2 ^(Pk + Ik) ~ E- The submanifold thus is 
a (2/ — 3)-dimensional sphere which we denote by 

4fflM = {(q- P) e 1^'^ : ^2(q, p) = , gi = PI = 0} . (6) 

This is a so called normally hyperbolic invariant manifold [9l[T0] (NHIM for short), i.e. 
■^NHIM invariant (since gi = j?i = implies qi = pi = 0) and the contraction and 
expansion rates for motions on S'j^him ^'^^ dominated by those components related 
to directions transverse to •S'^h^m- NHIM ^ can be considered to form the 

equator of the dividing surface (jS) in the sense that it divides into two hemispheres 
which topologically are (2/ — 2)-dimensional balls. All forward reactive trajectories 
(i.e. trajectories moving from reactants to products) cross one of these hemispheres, 
and all backward reactive trajectories (i.e. trajectories moving from products to 
reactants) cross the other of these hemispheres. Note that a trajectory is reactive 
only if it has Xi > (i.e. if it has sufficient energy in the first degree of freedom). 
Trajectories with Xi < are nonreactive, i.e. they stay on the side of reactants or on 
the side of products. 

Due to its normal hyperbolicity the NHIM has stable and unstable manifolds 
which are given by setting pi = —qi resp. pi = qi on the energy surface. Each of 
them have two branches. We denote the branches of the stable manifold by 

^^NHIM;r = {(q, P) ^ 1^'^ ■ ^^2(q, p) =E,pi = -qi > 0} , 
W^NHlM;p = {(q, P) e ■■ H2{CL, p) = E , pi = -qi < 0} , 



where is located on the reactants side of the dividing surface and is located 
on the product sides of the dividing surface. Similarly the unstable manifold has the 
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two branches 

W^NHIM;r = {(Q, P) ^ ■ -f^2(q, p) = ^, Pi = gi < 0} , 

W^NHlM;p = {(q, P) e E^'-^ : H2{q, p) =E,pi=qi>0}. 

The stable and unstable manifolds have the topology of spherical cylinders Rx S'^^^^. 
As they are of co-dimension 1 in the energy surface they can act as separatrices, 
dividing the energy surface into different components. In fact on the stable and 
unstable manifolds Ii is equal to zero. They thus lie 'between' the non-reactive 
(Xi < 0) and reactive (IiiO) trajectories. More precisely the forward reactive cylinder 

Cf = W^HIM;rUT^,:jHIM;p (9) 

encloses all forward reactive trajectories in the energy surface and separates them 
from all other trajectories. Similarly the backward reactive cylinder 

Cb = W^NHIM;pUW^^jHIM;r (10) 

encloses all backward reactive trajectories in the energy surface and separates them 
from all other trajectories. In this sense Cf and Cb form the phase space conduits of 
forward and backward reactions, respectively. 
Moreover, it is useful to define the lines 

DRPf = {(q,p) gR2/ . H2{q,p)=E,pi>0,pk = qk = 0,k = 2,...,f}, 
DRPb = {(q,p) G R'^ : i/2(q,p) = E, < 0, pk = Qk = 0, k = 2, . . . J} . 

These lines form the centerlines of the volumes enclosed by the forward and back- 
ward reactive cylinders Cf and Cb, respectively, and are therefore referred to as the 
forward and backward dynamical reaction paths, respectively. All forward reactive 
trajectories spiral about DRPf, and all backward reactive trajectories spiral about 
DRBb- Note that in the limit E ^ the forward and backward reactive cylinders Cf 
and Cb shrink to the forward and backward dynamical reaction paths, which in turn 
become the one-dimensional stable and unstable manifolds of the saddle equilibrium 
point in this limit. More precisely, if 

^saddle;r = {(q, P) G ^'^'^ ■ ^^2(q, p) = 0, Pf = -gi > 0, Pfc = % = 0, /c = 2, . . . , /} , 

^saddleip = {(q> P) G ■■ H2{ci, p) = 0, pi = -qi < 0, Pk = Qk = 0, k = 2, . . . , /} 

(12) 

denote the reactants and product branches of the stable manifold of the saddle, and 

W^saddle;r = {(q, P) ^ : i/2(q, p) = 0, pi = qi < 0, Pk = qk = 0, k = 2, . . . , f} , 

^saddle;p = {(q> P) ^ : -ff2(q, p) = 0, Pi = qi > 0, Pk = qk = 0, k = 2, . . . , f} 

(13) 



denote the reactants and product branches of the unstable manifold of the saddle 
then DRPf ^ I^|^,,,^, U W--,,,,^^ and DRPb ^ T^|^die;p U W^saddie;r for E ^ 0+. 
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For systems with / = 1 degree of freedom with a saddle, one can still separate 
the phase space (or energy surface) in a reactants and a products region in a similar 
manner as described above. However, most of the phase space structures defined 
above make no sense for / = 1. The case of one degree of freedom is special since in 
this case the trajectories are given by the level sets of the Hamiltonian. The question 
of whether a trajectory is reactive or not is thus completely determined by the (total) 
energy of the trajectory. 



2.2 The general (nonlinear) case 

For the general nonlinear case, consider a Hamiltonian function H which has an 
equilibrium point (a 'saddle') at which the corresponding linearized vector field has 



the same disposition of eigenvalues as in Sec. 2.1 In the neighbourhood of the saddle 



the dynamics is thus similar to that of the linear vector field described in Sec. 2.1 



In fact if follows from general principles that all the phase structures discussed in 



Sec. 2.1 persist in the neighbourhood of the saddle (which in particular implies that 
one has to restrict to energies close to the energy of the saddle). Moreover, these 
phase space structures can be constructed in an algorithmic fashion using a Poincare- 
Birkhoff normal form [6l [12]. Assuming that the eigenvalues w^, k = 2, . . . , /, are 
independent over the field of rational numbers (i.e. in the absence of resonances), 
the Poincare-Birkhoff normal form yields a symplectic transformation to new {normal 
form) coordinates such that the transformed Hamiltonian function truncated at order 
riQ of its Taylor expansion assumes the form 

Fnf(Xi,/2,...,//), (14) 

where Xi and Ik, k = 2, . . . , /, are constants of motions which (when expressed in 
terms of the normal form coordinates) have the same form as in ^, and H^F IS 
a polynomial of order no/2 in Ii and Ik, k = 2, . . . , f (note that only even orders 
riQ of a normal form make sense). The algorithm to compute this transformation is 
sketched in Appendix \K\ 

In terms of the normal form coordinates the phase space structures can be defined 
in a manner which is virtually identical to the linear case by replacing i?2(c[)P) by 



H^-pili, I2, . . . , //) in the definitions in Sec. 2.1 Using then the inverse of the normal 
form transformation allows one to construct the phase space structures in the original 
('physical') coordinates. The Poincare-Birkhoff normal form is therefore of crucial 
importance for the construction of the phase space structures governing reaction 
dynamics in a general Hamiltonian system with a saddle equilibrium point. 



3 Reduction of the N-body system 

We now want to study reaction type dynamics induced by the relative equilibria of a 
rotationally reduced A^-body system. We start by recalling the reduction procedure 
for the A-body system following [T7]. Let Xj, i = 1, A, be the position vectors of 
A^ bodies in . If a set of mass- weighted Jacobi vectors Sj,i = 1,...,A^— 1, are chosen 
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(for an example, see Sec. 5.2) and the center of mass of the system is assumed to be 
the origin, then the three degrees of freedom associated with overall translations are 
eliminated, and the kinetic energy takes a diagonal form 

1 ^ 

K = -Y^sl,. (15) 

i 

After a suitable choice of a body frame one can write Sj = R(0i, ^3)1"?, for some 
R'(^i) ^2, ^3) G SO{3), where 9i, 62, 63 is some set of Euler angles, and r^, i = I, ...,N— 
1, are the mass- weighted Jacobi vectors in the body frame. The can be expressed 
in 3A^ — 6 coordinates q^, called shape coordinates and their space is called the shape 
or internal space which we denote by M. Away from collinear configurations, as 
we will assume throughout, the translation reduced space is a fiber bundle over the 
shape space. For a conservative A/^-body system without external forces the potential 
energy can be viewed as a function on the shape space (i.e. it is a function of the 
shape space coordinates only). The goal now is to also express the kinetic energy (as 
far as possible) as a function of the shape space coordinates. To this end we make 
the following definitions. 

Let R G SO{3) denote the rotation from the center of mass frame to the body 
frame and L denote the angular momentum 

N 

L = ^SiXSi. (16) 

i=l 

Then the body velocities and body angular momentum are defined, respectively, by 

h = R^s„ (17) 

and 

J = R^L . (18) 
The moment of inertia tensor M{q) has the components 

iV-l 

Mijiq) = J2 ("^^"^'J ~ '^kinj) , (19) 
k=l 

where = (r^i, rk2,fkz) in body coordinates, and is the so called gauge potential 

A,{q)=M-\q).Y^^i^r,x^J . (20) 
After defining the metric 

g,. = P^^-A,.M.A, (21) 
dqf, dqu 

(where we here and in the following use the Einstein convention of summation over 
repeated indices which run from 1 to 3A^ — 6) the kinetic energy becomes 

K = ^(^ + ^M^At) - M- {uj + Ayqy) + ]^g^yq^qu . 



Here uj is the angular velocity which is the vector corresponding to the skew-symmetric 
matrix i^'^i? by the natural isomorphism 



(22) 
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By using the equation 



dK 

i=-^=M{u + A^q^), (23) 



the conjugate momenta of the shape space coordinates are obtained by 

dK 

Pf, = -qT- = 9fiuqv + J • . (24) 
The fully reduced Hamiltonian is then given by 

= • • i+^g^'^'ip^ - J • A^){p, - J • A,) + V{qu ...,q3N-e) • (25) 

Here the first term is called rotational or centrifugal energy and the second one is 
called vibrational kinetic energy. 

The equations of motion are obtained to be 

OH 

q^ = dH/dp^, Pf, = -dH/dqf,, j = Jx— (26) 



for /U = 1, ...,3A^ — 6, [12]. The last equation is equivalent to [H] 

ja = {Ja,H}, a =1,2,3, (27) 

where J = (Ji, J2, J3). 

The theory underlying the reduction described above is the orbit reduction method 
of Marie [33]. Let Q denote the translational reduced configuration space of non- 
collinear configurations. Then the action of 50(3) on that space is free which implies 
that the lifted action of S0{3) on T*Q is also free besides being symplectic. The 
angular momentum can be considered as a map (a so called momentum map) J from 
T*Q to the dual space of the Lie algebra of S0{3), 50(3)*, which is identified with 
R^. After writing the Hamiltonian in terms of a body fixed frame, it takes an 50(3)- 



invariant form as in (25). As the (space fixed) angular momentum is constant, say 
and the magnitute of the body angular momentum is also a constant of motion, the 
dynamics can be reduced to the inverse image under the angular momentum map of 



the body angular momentum sphere Or = Su^u which is the coadjoint orbit of the 



50 (3) action. Then (25) may be considered as a function on 7^^(0^)750(3) which 
is symplectomorphic to the Marsden-Weinstein space J~^(C)/50(2), where 50(2) is 
the isotropy subgroup of the coadjoint action of 50(3) on 50(3)*. 
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4 Canonical coordinates on the reduced phase 
space 



It can be seen by (25) that {q^,Pf_i,Ja) are some coordinates on the reduced space 
where n = 1, .., 3A^ — 6 and Ja,a = 1, 2, 3, are the components of the body angular 
momentum vector J. Note here that {q^,p^) are canonical variables whereas the J a 
are not canonical. If one uses coordinates = — J • in place of p^, then 
il^iiVfi) become non-canonical too [T7] but this way it easy to see that the reduced 
space is locally diffeomorphic to the product of the angular momentum sphere and 
the cotangent bundle of the internal space: S'pn x T*M. This is in fact a general 
fact for cotangent bundle reduction, known as fibration cotangent bundle reduction 
[HJ. It is to be noted here that the symplectic structure on the reduced space is 
the sum of ±1/ || J|| times the volume form of the sphere, the canonical symplectic 
structure on T*M and some magnetic terms 03]. The magnetic terms come from 
the gauge potential, and are directly related to the non-Euclidean structure of the 
shape space. The goal now is to use some canonical coordinates u, v on Spy instead 
of the coordinates Ja- The related symplectic structures on the coadjoint orbit Sr, 
where r = ||J|| , are 

w± = ±-dA (28) 
r 

where dA is the volume form on the sphere, or in any coordinates V on 

w± = ±--\/det \guv]du A dv , (29) 
r 

where [guv] is the matrix of the induced Riemannian metric on S'^ by the Euclidean 
metric on R^. The canonical coordinates with respect to this symplectic structure 
are those that give 

w± = ±du A dv . 
One possible choice to achieve this is [l5] 

Ji = \/ r'^ — cos n, J2 = \/ r"^ — v"^ sinn, J3 = v , (30) 

where < u < 2tt, —r < v < r (see Fig. [Tj^a)). Let u, v be some canonical coordinates 
for the positive signed symplectic structure w+ and 1 < fi < 3N — 6, be any 

choice of shape coordinates and their conjugate momenta, respectively. In terms of 
these coordinates, after relabeling u = qo,v = po, and setting z = (zj) = {qi,Pi),i = 
0, 3A'' — 6, the equations of motion read 

ii = {z„H} (31) 

or equivalently 

qi = dH/dpu Pi = -dH/dqi (32) 

for i = 0, 3A^ — 6, where the Poisson bracket on the reduced space has the standard 
form 

If \ = S^3N -6—— - —— (33) 
^ dqi dpi dpi dqi ' 
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(b) 




Figure 1: {u,v) coordinate lines on the angular momentum sphere for (a) the definition 



of {u,v) in pOj) and (b) the definition of {u,v) in (38). 



This canonical form of the equations of motion has many advantages. Although 
the coordinates p^, fj, = 1, 3A^ — 6, are gauge dependent, for examining the local 
behaviour of the reduced system they can be used to study, for example, the stability 
of relative equilibria \4:1\ I42j . We will use the canonical equations of motion to 
construct the phase space structures governing reaction type dynamics associated 
with saddle type relative equilbria in the reduced system from a Poincare-Birkhoff 
normal form in the same way as described in Sec. [2} 



Relative equilibria are trivial solutions of the equations of motion ( 26 ) or equiv^ 



reduced equations 



alently (|32|). Using (|26|) one finds that the relative equilibria are the solutions of the 

J X (M^^ - J) = 0, (34) 

P^^ = J-A^, (35) 

d 

(J-M-i • J+y(g)) = 0. (36) 



dq. 



This implies that at relative equilibria the body angular momentum J is parallel to 
a principal axis, i.e. the molecule is rotating about one of its principal axis. Using 
this fact the relative equilibria can be found from the critical points of the effective 
potential 

Ves{q)3 ■ • 3+V{q) , (37) 

where J is a fixed vector of a given modulus r parallel to a chosen principal axis. Some 
of these relative equilibria are related to the equilibria of the system with zero angular 
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momentum. In fact for a triatomic molecule without any discrete symmetries, six 
families of relative equilibria are born out of any generic equilibrium if the modulus of 
the angular momentum is increased from zero. Here 'generic' means that the rotations 
act on the equilibrium configuration freely (this excludes collinear configurations) 
[42j . We will follow this procedure to find relative equilibria in the example given in 
Sec. EH 

Given a relative equilibrium we can determine its (linear) stability using the equa- 



tions of motion (32 ). If the stability is of saddle type (as defined in the introduction) 



we can apply the Poincare-Birkhoff normal form procedure described in Appendix [A| 



to the equations of motion ( 32 ) in order to construct the phase structures governing 
reaction type dynamics as described in Sec. [2j To this end we note that {u, v) defined 
as in (30) are singular at the north pole and south pole J = (0, 0, ibr) of the angular 
momentum sphere. The coordinates are therefore not useful to study relative equi- 
libria associated with rotations about the third principal axis [l6] . In order to study 



these relative equilibria one can redefine {u, v) according to 

Ji = V, J2 = -y/r^ — ti2 sinn, J3 = \/ — cos u (38) 

in which case the coordinate singularities are at the points J = (0, ±r, 0) on the 
angular momentum sphere (see Fig. [T|b)). 

5 Examples 

In the following we illustrate the approach above for the limiting case of a rigid 
molecule and the case of a triatomic molecule. 

5.1 The limiting case of a rigid molecule 

In the limiting case of a rigid molecule the internal degrees of freedom are frozen, i.e. 



Qfi = Qfj. = 0, 1 < fj. < 3N — 6. In this case, one can neglect the potential in (25) and 
the Hamiltonian reduces to 

H =^3 ■ M"^ • J (39) 

which describes the free rotation of a rigid body. Choosing the body frame to coincide 
with the principal axes, M (and its inverse M~^) becomes diagonal and (39) becomes 



1 t2 t2 t2 

where Mi, M2 and M3 are the principal moments of inertia. In this case the angular 



momentum part of (26) become the classical Euler equations [15]. In terms of the 



coordinates {u,v) defined in (30) these equations become 



u = dH/dv, V = -dH/du . (41) 

Here the role of v is seen to be a momentum conjugate to v but when we choose the 
other symplectic form u}~ = —-dA (see (28)), then the roles of u and v are reversed. 
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If we choose the coordinates ( 30 ) , then the Hamiltonian becomes 



1 / (r^ — t!^) cos u (r^—v)smu v \ 

^(^'^)= — d — ^ — d — + 1^) 



So, by (41) the equations of motion are obtained to be 



V cos^ u V sin^ u v 
+ 



Ml Ms Ms ' 

(r^ — f^) sinncosn (r^ — f-^) sin?xcos?x 
^ ~ Ml M2 ' 

In the case of the Euler top where the moments of inertia are mutually different, i.e. 
where we can assume without restriction that Mi < M2 < M3, there are six relative 
equilibria given by the points (itr, 0, 0), (0, ±r, 0) and (0, 0, ±r) on the J sphere. The 
dynamics near the first two can be studied in terms of the coordinates {u, v) defined 
according to ( pO| ) where they correspond to the points (0, 0) (for J = (r, 0, 0)), (tt, 0) 
(for J = (-r,0,0)), (7r/2,0) (for J = (0,r,0)) and (37r/2,0) (for J = (0,-r,0)) 
(see Fig. [2]^a)). The dynamics near the relative equilibria J = (0,0, itr) (and again 



J = (0, zbr, 0)) can be studied in terms of (n, v) defined according to (30) where they 
correspond to the points the points (0, 0) (for J = (0, 0, r)), (vr, 0) (for J = (0, 0, — r)), 
(7r/2,0) (for J = (0,r,0)) and (37r/2,0) (for J = (0,-r,0)) (see Fig.[2][b)). The re- 
duced Hamiltonian has local maxima at J = (itr, 0, 0) which correspond to rotations 
in either direction about the principal axis with the smallest moment of inertia, and 
minima at J = (0, 0, ibr) which correspond to rotations in either direction about the 
principal axis with the largest moment of inertia. The rotations about the principal 
axis with middle moment of inertia corresponding to J = (0, itr, 0) are unstable. The 
corresponding relative equilibria are of saddle type. The eigenvalues associated with 
the linearized vector field are ±y^MiM3(M2 - Mi)(M3 - M2)/(MiM2M3). The 
eigenvectors are shown in Fig. [2] Since the reduced rigid body motion has one 
degree of freedom the reaction dynamics associated with these saddles is trivial (see 



the remarks at the end of Sec. 2.1). One can divide the phase space of the reduced 



system into the regions v < and v > in terms of the coordinates (30) or equiva- 
lently in the regions |m| < tt/2 and \u\ > 7r/2 in terms of the coordinates (38) (where 
we used the periodicity of u in the last case; see Fig.[2|. For energies below rV(2M2) 
which is the energy of the saddle one cannot reach one region from the other. For 
energies above r^/(2M2) however this becomes possible. 

5.2 Triatomic molecules: HCN 

For a second example, we consider a triatomic molecule. We will use the xxy gauge 
and Dragt's coordinates as follows [T7j. Let xi,X3,X3 be the position vectors of 3 
bodies in R^. If mass- weighted Jacobi vectors si,si are chosen as 

/—f \ mixi + 7713x3 

Sl = Vm(xi - X3) , S2 = VM2(X2 , ), (42) 

ra\ + m3 

where 

mim3 m2imi+'m3) 

fii = and = (43) 

nil + m3 mi + m2 + m^ 
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(b) 



V 





Figure 2: Contours of the reduced rigid body Hamiltonian in the {u, v) coordinate plane 
with (m, f) defined according to (30) (a) and with {u,v) defined (38) (b). The energy- 
increases from hght to dark shading. The bold lines and arrows indicate the eigenvectors 
and the corresponding directions of the Hamiltonian fiow for the saddles. The parameters 
are Mi = 1, Ma = 2, M3 = 3 and r = 1. 




Figure 3: Definition of the Jacobi vectors Si and S2 and the corresponding angle (p (see 
(42) and (5.2)) for the HCN molecule. Carbon is associated with the first mass, nitrogen 
is associated with the second mass and hydrogen is associated with the third mass. 
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are the reduced masses then the corresponding Jacobi coordinates {pi , P2,(p) are 
defined by 



Pi 



|si| 



P2 



Si • S2 = P1P2 cos I 



where < (p < it. From the Jacobi coordinates one can define the coordinates 

{wi,W2,W3) = {pj - /)2,2/OlP2 COS 0,2/91/92 sine/)) , (44) 

where 'Wi,W2 G K and W3 > 0. Dragt's coordinates {uj,x,'4') are now polar coordi- 
nates in the {wi, 11)2,103) coordinate space: 



{■wi,W2,W3) = {uj COS X COS ^p,UJ COS X simp, CO sin x) , 



(45) 



where w > 0, < x ^ 7r/2, < -0 < 2tt. Note that ^ is the latitude, not the colatitude. 
The inertia and metric tensors are diagonal in Dragt's coordinates |17j : 



M 



cj sm 





2 X 



uj cos^ I 



and 



4uj 
^ 

UJ 





w cos^ X 



The gauge potential becomes 
A, 



= (0,0,0), Av = (0,0,--sinx). 



Hence, if the coordinates (30) are used then after relabeling u 
QiiX = Q2,ip = qs, the reduced Hamiltonian becomes 



qo,v 



(46) 
(47) 
(48) 



H{q,p) 



(r^ - Pq) cos^ ( 
2qi sin2 f 

+^(91,92,93) 



+ 



(r^ — Pq) sin^ I 
2qi cos2 f 



i -Po , 9„ ^2 , (2p3+po smg2) 

+ o ^ ^^l^'l ^ ^ 2 

2gi gi 2^1 cos^ q2 



or if the coordinates ( 38 ) are used one has 

[r -Po)sm go 



H{q,p) 



+ 



2qi cos2 f 



2gi sin^ f 

(2p3 + - p2 cos qo sin 92)^ , ^ 

+ S + ^(^1' 92, 93 

2gi cos^ q2 



2gi qi 



As a concrete example we consider the HCN molecule where we use the potential 
by Murrell, Carter and Halonen |47j . The molecule exists in the form of two isomers 
which correspond to the coUinear configurations HCN (hydrogen cyanide; the corre- 
sponding equilibrium has energy -13.5914 eV) and CNH (hydrogen isocyanide; the 
corresponding equilibrium has energy -13.1065 eV). For zero angular momentum, the 
reaction dynamics of the isomerization from HCN to CNH or vice versa is induced by 
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Figure 4: Potential energy lines (grey) of the Murrell, Carter and Halonen HCN 
potential in the section dcN = 1.1394 A together with (a) the projection of the real 
eigenvectors of the saddle equilibrium of the system with vanishing angular momentum 
and (b) the path traced by the relative equilibrium under variation of r projected to this 
plane for rotations about the first (full line), second (dashed line) and third (dotted line) 
principal axis. 



a saddle equilibrium point which corresponds to a triangular (noncoUinear) configu- 
ration (and has an energy of -12.0827 eV). The internuclear distances corresponding 
to this saddle are 

dcN = 1.1394 A , dcR = 1.1206 A , c^nh = 1.4950 A . (49) 

The eigenvalues of the matrix associated with the linearized vector field at this equi- 
librium are ±589.311, ±418.691 and ±213.08 (in units of inverse pico seconds) indi- 
cating that it is of saddle x center x center stability type. 

We assign the Dragt's coordinates by identifying the first body with carbon, the 
second body with nitrogen and the third body with hydrogen (see Fig. [s]). The 
internuclear distances are then given by 

dcN = i — Vqi + Qi cosg2Cosg3, 



^CH = il :^{qi + qicosq2COsq3) + ^{qi - qicosq2COsq3) '^^ qi cos g2 sings, 

2mf 2/^2 rniy/JPz 



dNH = \l :p^{qi + qicosq2COsq3) + ^{qi - qicosq2COsq3) -\ ^^gi cos 92 sings 

V 2m^ 2/X2 "^3^/Ai2 

where mi, m2 and ms are the masses of carbon, nitrogen and hydrogen, respectively. 



and m and ^2 are the reduced masses defined according to (43). The internuclear 
distances can also be viewed a coordinates on the shape space. The equipotential 
lines of the Murrell-Carter-Halonen potential intersected with the plane dcN = 1-1394 
(which contains the saddle) are shown in Fig. |4]^a). In Fig. Qa) also the projection 
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r 



Figure 5: Energies of the relative equilibria as a function of the modulus of the angular 
momentum r. The energy is measured in eV and r is measured in units of h. For fixed r, 
the relative equilibrium involving rotations about the principal axis with smallest (highest) 
moment of inertia has highest (smallest) energy. 



of the eigenvectors corresponding to the real eigenvalues of the saddle are shown. 
Due to the time reversal symmetry of the system without angular momentum they 
project to the same line. 

For the zero angular momentum saddle, the moment of inertia tensor has the 
diagonal elements [1.1645, 8.6470, 9.8115] (in the units of xeV). Due to the similar 
masses of carbon and nitrogen the values of the second and third moment of inertia 
are quite similar. We consider the six relative equilibria resulting from this saddle as 
the modulus r of the angular momentum is increased from zero (see Sec. [4]). Since 
the directions of rotations about either principal axis are related by symmetry there 
are effectively three families. We numerically compute these relative equilibria as a 
function of r from the critical points of the effective potential Veff((7i, (/2; 93) in (37) 
where the angular momentum J is (r, 0,0), (0,r, 0) or (0, 0,r), respectively, using a 
Newton procedure. We show the energies of these relative equilibria as a function of 
r in Fig. [5| Since for a rigid body, the rotation about the principal axis with middle 
moment of inertia is unstable it is to be expected, that at least for small values 
of r, the two of the relative equilibria involving rotations about the first and the 
third principal axis are saddle x center x center x centers and the relative equilibrium 
involving rotations about the second principal axis is a saddle x saddle x center x center 
(which has two pairs of real eigenvalues and two complex conjugate pairs of imaginary 
eigenvalues). Linearizing the vector fields at the relative equilibria we find that this 
is indeed the case. In Fig. [6]we show the disposition of the eigenvalues in the complex 
plane as a function of r for the three cases. For r — )• 0, the eigenvalues reduce in each 
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case to the eigenvalues of the equihbrium in the system without angular momentum 
plus a double eigenvalue at zero. 

Figure |4] shows the projections of the paths traced by the relative equilibria under 
variation of r to the {den, c^nh) plane. One finds that in this projection the relative 
equilibrium corresponding to rotations about the first principal axis is moving much 
faster than the other two relative equilibria as r is varied. Note that the position 
of the relative equilibria on the angular momentum sphere does not change under 
variation of r. 

The saddle type relative equilibria corresponding to rotations about the first and 
the third principal axis induce reaction type dynamics. The dynamics and its physical 
implications will be discussed in another publications. We here restrict ourselves to 
a brief discussion of the reaction paths in the limit where the energy approaches 
the energy of the relative equilibrium from above (see the discussion in Sec. [2]). 
In this limit the reaction paths reduce to the one-dimensional stable and unstable 
manifolds of the relative equilibria. At the relative equilibria these lines are tangent 
to the eigenvectors corresponding to the real eigenvalues of the linearization of the 
reduced system. We therefore show pictures of the pairs of real eigenvectors of the 
relative equilibria in Fig. [7| For the rotation about the first principal axis, the 
projections of the eigenvectors corresponding to positive and negative eigenvalues 
to the (dcHi'^NH) plane are almost identical (similar to the angular momentum zero 
case). Their projections to the {u, v) plane however, are quite different (see Figs. [Tj^a) 
and (b)). As the example indicates the relative equilibrium (in its projection to 
the shape space) might be quite far away from the position of the corresponding 
equilibrium of the system without angular momentum, and hence the reaction paths 
change quite considerably as a function of r. In the case of rotations about the third 
principal axis the projections of the eigenvectors of the two real eigenvalues to the 
{dcH,djqii) plane are quite different whereas their projections to the {u,v) plane is 
zero. The reason for the latter is related to the fact that in the case of rotation about 
the third principal axis the motion is planar (in the xxy gauge we are using the 
angular momentum vector is perpendicular to the plane spanned by the three bodies 
for rotations about the third principal axis). Although this relative equilibrium 
(in its projection to the shape space) is still relatively close to the corresponding 
equilibrium of the system without angular momentum the different directions of the 
eigenvectors indicate that the forward and backward reaction paths might globally 
be quite different for the system with and without angular momentum. 

6 Conclusions and Outlook 

In this paper we gave an explicit formalism to locally construct canonical coordinates 
in rotational symmetry reduced A^-body systems. This allows one to study relative 
equilibria in such systems and use Poincare-Birkhoff normal form to determine the 
phase space structures that govern reaction type dynamics near relative equilibria of 
saddle stability type. We briefly illustrated the formalism by applying it to the saddle 
relative equilibria which induce the isomerization reaction dynamics of HCN/CNH 
with rotation. We here restricted ourselves to the study of the reaction paths near 
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Figure 6: Paths traced by the eigenvalues (in units of inverse pico seconds) of the hnearized 
vector fields as the modulus of the angular momentum r is varied for rotations about the 
first principal axis (left panel), the second principal axis (middle panel) and the third 
principal axis (right panel). The ranges for r are about the same as in Fig. [s] 

the saddle relative equilibria. A more detailed study of the isomerization dynamics 
which includes a Poincare-Birkhoff normal computation will be presented elsewhere. 
As an outlook we give a short list of open questions and future directions related to 
the work presented in this paper. A question that we did not discuss is the recon- 
struction problem, i.e. the study of what the corresponding motions of the reduced 
system are in the full original system. Especially from the perspective of reaction 
dynamics this is an interesting and important problem. The formalism introduced 
in this paper is not restricted to the study of relative equilibria of saddle type. Re- 
cently also the relevance of higher rank saddles for reaction type dynamics has been 
studied |48| HH] • In the HCN / CNH isomerization problem with angular momentum 
a saddle X saddle X center X center relative equilibrium (and its symmetric partner) is 
presented by the rotation about the principal axis with the middle moment of inertia 
(and the two directions of rotation). Interestingly, for a fixed value of the modulus 
of the angular momentum, these rank 2 saddle relative equilibria are quite close in 
energy to the saddle relative equilibria corresponding to motion about the principal 
axis with the largest moment of inertia (see Fig. [s]) . It would be interesting to study 
the interplay between the phase space structures associated with the rank two saddles 
and usual saddles and their influence on the global dynamics. In [501 ED IS2] (see also 
[53]) a procedure based on transition state theory to determine the microcanonical 
volume of initial conditions that lead to reactive trajectories has been introduced. 
The formalism presented in this paper allows one to study such reactive volumes 
for rotational symmetry reduced systems. An interesting question for applications 
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(a) (b) 




Figure 7: Saddle type relative equilibria and the corresponding real eigenvectors projected 
to the {u, v) plane and {den, c^nh) plane, respectively for rotations about the first principal 
axis with r = 33.5 , /i ((a) and (b)), and for rotations about the third principal axis with 
r = 50,n ((c) and (d)). 
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would be to study how this reactive volume changes when the modulus of the angular 
momentum is varied. The formalism presented in this paper at first only applies to 
relative equilibria of noncollinear configurations. We remark however that it can be 
quite simply adapted to also include collinear configurations. For this purpose, one 
can consider the Hamiltonian and coordinates in [l2] and choose suitable coordinates 
in place of the angular momentum. In this paper we did not discuss the quantum 
reaction dynamics of rotating molecules. In |54^ [T2] the Weyl symbol calculus has 
been used to develop a quantum mechanical analogue of the Poincare-Birkhoff normal 
form which provides an efficient procedure to compute quantum reaction rates and 
the related Gamov-Siegert resonances. The same procedure can in principle also be 
applied to the relative equilibria of rotational symmetry reduced molecules. However, 
the reduction of the quantum A^-body system has some differences to the classical 
case. For example, for the quantum three body problem one can separate rotations 
and vibrations whereas this is not possible in the classical case [55] . It would be very 
interesting to study the quantum reduction problem from the perspective of reaction 
dynamics. Finally we mention that shortly before the submission of our paper a re- 
lated paper on the phase space structures governing reactions in rotating molecules 
was published [56J- That approach is based on a different normalization procedure 
which is treating the body angular momentum in a non-canonical manner using a 
rigid body approximation. It would be interesting to compare that approximation 
with the study in this paper in more detail. 



Acknowledgments 

V.Q. is grateful for the hospitality of the Johann Bernoulli Institute for Mathematics 
and Computer Science at the University of Groningen where the research for this 
paper was carried out, and to Namik Kemal University for financial support during 
his stay. 



Appendix 

A The Normal Form Algorithm 

We here describe the algorithm to compute the Poincare-Birkhoff normal form in the 



neighbourhood of a saddle equilibrium point of the form described in Sec. 2.2 The 
implementation of this algorithm in a computer program is freely available [57] . 

Let H2 denote the quadratic Hamiltonian which gives the linearized Hamiltonian 
vector field at the saddle. One says that a Hamiltonian H is in normal form if H 
Poisson commutes with its quadratic part, i.e. 

fc=l 



In general H is not in normal form. However, for any given order no of the Taylor 
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expansion of H one can find a symplectic transformation to new phase space coordi- 
nates in terms of which the transformed H truncated at order uq is in normal form. 
This symplectic transformation is constructed from a sequence of the form 

{qupi) = z = z(°) ^ z(i) ^ z(2) ^ z(3) ^ . . . ^ z("«) , (51) 

where z^"^ is obtained from z^*^"^-* by means of a symplectic transformation 

^ = </,v^„z('^-^) . (52) 

generated by a polynomial Wn{z) of order n, i.e. 

. 1^,1 ^ifli _„| ^ (53) 



Wn e := span jg^^i . . . ...py : \a\ + 



Here |q| = I]fc=i 



Y^i=il3k- More precisely, the in (52) denote the 



time-one maps of the flows generated by the Hamiltonian vector fields corresponding 
to the polynomials Wn (see [12] for the details). The maximum order no in (51 ) is the 
desired order of accuracy at which the expansion will be terminated and truncated. 

Expressing the Hamiltonian H in the coordinates z^"'\ n = 1, . . . ,no, we get a 
sequence of Hamiltonians H^^^ , 

H = ^ ^ ^ ^(3) ^ _ _ , ^ ^(no) ^ (54) 
where forn = 1, . . . ,no, if(")(z(")) = f("-i)(z("-i)) = ij("-i)((/>^^^z(")), i.e. 

^(n)^^(n-l)^^-l^_ (55) 

To avoid a proliferation of notation we will in the following neglect the superscripts 
(n) for the phase space coordinates. 

In the first transformation in (51) we shift the equilibrium point zg to the origin, 
i.e. z I—)- := z — zq. This gives 

i7(i)(z) = i/(o)(z + zo). (56) 

The next steps of the normal form procedure rely on the power series expansions of 



(57) 



s=2 



(n) 

where the Hs are homogenous polynomials in W": 



E 



H 



(n) 

aiv.a/,/3iv,/3/ 

ai! . . . a/!/3i! . . . 



■9/ Pi •••P/ 



(58) 



For n = 1, the coefficients in (58) are given by the Taylor expansion of about 
the origin. 

' z=0 



H 



(1) 
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For n > 3, the coefficients in (58) are obtained recursively. For n = 2, i.e. the second 
step in the sequence of transformations (51), the coefficients in (58) are determined 
by a Unear transformation of the phase space coordinates according to 

z ^ 4>W2{z) := M z . (60) 

Here, M is a symplectic 2/ x 2/ matrix which is chosen in such a way that the 
transformed Hamiltonian function 



H^'^\z) = h'^^^M-^z) 



(61) 



assumes the same form as in ([T]). Section 2.3 of Ref. [12j provides an exphcit procedure 
for constructing the transformation matrix M. 

For the first two steps in the sequence (51), we actually did not give explicit 
expressions for the generating functions Wi and W2- For conceptual reasons (and to 
justify the notation) it is worth mentioning that such expression can be determined 
(see P2j)- The next steps in (51) though rely on the explicit computation of the 
generating functions Wn with n > 3. To this end it is convenient to introduce the 
adjoint operator associated with a phase space function A: 



adA-. adAB = {A, B} . 



(62) 



The transformation (52) then leads to a transformation of the Hamilton function 
^(n-i) with n > 3 which in terms of the adjoint operator reads 

oo ^ 

^^"^=Efc![^dH.j'^("-^^. (63) 



fe=0 



In terms of the Taylor expansion defined in Eqs. (57 59) the transformation intro- 
duced by Eq. ( 63 ) reads 



Ln-2j ^ 



r{"-l) 

s-fc(n-2) ' 



(64) 



k=0 



where [-J gives the integer part of a number, i.e., the 'floor '-function. 



Using Eq. (64) one ffiids that the transformation defined by (63) satisfies the 
following important properties for n > 3. Firstly, at step n, n > 3, the terms of order 
less than n in the power series of the Hamiltonian are unchanged, i.e. 



= iJ^-i) , for s<n. 



so that, in particular, H2^^ = H"!^' . Defining 



r{2) 



'D = ad 



H. 



(2) 



(65) 



(66) 



we get for the term of order n. 



(67) 
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This is the so-called homological equation which will determine the generating func- 

(n) (n) 

tions Wn for n > 3 from requiring "DHn = 0, or equivalently Hn to be in the 



kernel of the restriction of V to W". In view of (67) this condition yields 



ffi^n^'^) -VWn G KerP| 



Section 3.4.1 of Ref. [12] provides the explicit procedure of finding the solution of 
Eq. ( 68 ) . In the generic situation where the linear frequencies a;2 , . . . , in ([T| are 



rationally independent, i.e. m20J2 + ■ ■ ■ + rufCOf = implies m2 = . . . = m/ = for 

(n) 

all integers m2, . . . ,mf, it follows that for odd n, Hn = 0, and for even n, 

e span |x"i/2"'^3"' ■■■If^'-\a\=n/2^, (69) 
where I = {pj - ql)/2 and h = {ql + p|)/2, with k = 2,...,f. 



Applying the transformation (63), with the generating function defined by (67), 



for n = 3, ... ,77-0, and truncating the resulting power series at order no one arrives 
at the Hamilton: 
Hamiltonian H: 



at the Hamiltonian H^p^ corresponding to the tiq^ order normal form (NF) of the 



<j?)(z) = i?o + 5^//i"")(z). (70) 



s=2 

We stress that H^p '^ represents an order approximation of the original Hamil- 
tonian H obtained from expressing H in terms of the normal form coordinates znf 
which in turn are obtained from the symplectic transformation of the original coor- 
dinates z = {qi,Pi) 



ZNF = 4>{'^) = i't'Wno ° '/'Wno-l O • • • O 4'W2 ° (f^Wl) i^) ■ (71) 

This is why it is legitimate to use H^p^ instead of H in analyzing the dynamics in 
the neighbourhood of the saddle. 

We emphasize that the full procedure to compute H^p^ and the corresponding 
coordinate transformation is algebraic in nature, and can be implemented on a com- 
puter. A computer program is freely available from [57 \ . 
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